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Heavy quark and charmonium production as well as their space-time evolution are studied in 
transport simulations of heavy-ion collisions at RHIC and LHC. In the partonic transport model 
Boltzmann Approach of MultiParton Scatterings (BAMPS) heavy quarks can be produced in ini- 
tial hard parton scatterings or during the evolution of the quark-gluon plasma. Subsequently, they 
interact with the medium via binary scatterings with a running coupling and a more precise Debye 
screening which is derived from hard thermal loop calculations, participate in the flow and lose 
energy. We present results of the elliptic flow and nuclear modification factor of heavy quarks and 
compare them to available data. Furthermore, preliminary results on suppression at forward 
and mid-rapidity are reported for central and non-central collisions at RHIC. For this, we study 
cold nuclear matter effects and the dissociation as well as regeneration of 7/ in the quark-gluon 
plasma. 



XLIX International Winter Meeting on Nuclear Physics 
24-28 January 2011 
BORMIO, Italy 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Heavy quarks and charmonium at RHIC and LHC 



Jan Uphoff 



1. Introduction 

Several experimental observations indicate that a deconfined medium consisting of quarks and 
gluons - the quark-gluon plasma (QGP) - is produced in ultra-relativistic heavy-ion collisions 
[1, 2]. Charm and bottom quarks are an ideal probe for the early stage of these collisions since 
they can only be created in initial hard parton scatterings of nucleon-nucleon interactions or in the 
medium in the beginning of the QGP phase, where the energy density is still large. After then- 
production they interact with other particles of the medium and can, therefore, reveal important 
information about the properties of the QGP. Flavor conservation renders them as a unique probe 
since they are tagged by their flavor even after hadronization. 

The experimental data of the elliptic flow V2 and the nuclear modification factor Raa of heavy 
quarks [3-5] show that the energy loss of charm and bottom quarks is comparable to that of light 
quarks. Whether this large energy loss is due to collisional or radiative interactions - or both (or 
even other effects) - is under investigation (see [5] for a recent overview and comparison with 
data). 

J/Y suppression in heavy-ion collisions was proposed to be a signature of the QGP a long 
time ago [6], but it remains challenging to disentangle the contributions of hot and cold nuclear 
matter effects to the measured suppression (see e.g. Ref. [7] for a recent overview). Since lattice 
calculation indicate that J/y can survive in the QGP to some extent [8, 9], a partonic transport 
model such as BAMPS is an ideal framework to investigate the influence of J/y dissociation and 
regeneration on the J/y suppression. 

This article is organized as follows. After the introduction of the parton cascade BAMPS we 
will discuss the production of heavy quarks at RHIC and LHC. In Sec. ^ our results on the elliptic 
flow and nuclear modification factor of heavy quarks at RHIC are compared to the experimental 
data. Furthermore, preliminary results on J/y suppression are presented in Sec. || and in Sec. ^ we 
conclude with a short summary. 



2. Parton cascade BAMPS 

For the simulation of the QGP phase we use the partonic transport model BAMPS [10, 11], 
which stands for Boltzmann Approach of MultiParton Scatterings. BAMPS simulates the fully 
3 + 1 space-time evolution of the QGP produced in heavy ion collisions by solving the Boltzmann 
equation. 



dynamically for on-shell partons with a stochastic transport algorithm and pQCD interactions. % 
are the relevant collision integrals, and /i(r,p,,f) the one-particle distribution function of species 
/ = g, c, c, b, b,J/Y, since light quarks are not included yet. In addition to the binary collisions 
2 — )• 2, also 2 -f-)- 3 scatterings for the gluons are possible. That is, the following processes are 
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Details of the model and the employed cross sections can be found in [10-12]. 

3. Heavy quark production at RHIC and LHC 

In heavy-ion collisions charm and bottom quarks are produced in hard parton scatterings of 
primary nucleon-nucleon collisions or in the QGP. To estimate the initial heavy quark yield, we 
use PYTHIA [13] and scale from proton-proton collisions to heavy-ion collisions with the number 
of binary collisions. Secondary heavy quark production in the QGP is simulated with BAMPS. 
For the initial gluon distributions, the mini-jet model, the color glass condensate model and also 
PYTHIA are employed. 

In Au-i-Au collisions at RHIC with ^snn = 200 GeV between 0.3 and 3.4 charm pairs are 
produced in the QGP, depending on the model of the initial gluon distribution, the charm mass and 
whether a K = 2 factor for higher order corrections of the cross section is employed [12]. This 
is only a small fraction of the initially produced charm quarks and can be neglected for the most 
probable scenarios. 

At LHC with the much larger initial energy density, secondary charm production is enhanced 
and not negligible. In Pb-i-Pb collisions with ^snn = 5.5 TeV between 1 1 and 55 charm pairs are 
produced during the evolution of the QGP [12] (see left panel of Fig. |l|). These values are of the 
same order as the initial yield. As is shown in the right panel of Fig. even in the 2010 run with 
y/si^ = 2.76 TeV between 5 and 28 charm pairs are created during the QGP phase. 

Bottom production in the QGP, however, is very small both at RHIC and LHC and can be 
safely neglected. As a consequence, all bottom quarks at these colliders are produced in initial 
hard parton scatterings. 

Further details on heavy quark production can be found in Ref. [12]. 



4. Elliptic flow and nuclear modification factor of heavy quarks at RHIC 

The elliptic flow and the nuclear modification factor 

l pl-Py \ „ SNaa/AptAj 
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Figure 1 : Number of charm quark pairs produced in a central Pb+Pb colli sion at LHC with y^s^N = 5.5 TeV 
(top) and y^SNN = 2.76TeV (bottom) according to BAMPS. The initial parton distributions are obtained with 
PYTHIA and the mini -jet model (and the color glass condensate for ^snn = 5.5 TeV). In all cases the initial 
charm quarks are sampled with PYTHIA for better comparison. 
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(Px and py are the momenta in x and y direction in respect to the reaction plane) of heavy quarks 
at mid-rapidity are observables which are experimentally measurable and reflect the coupling of 
heavy quarks to the medium. A large elliptic flow and a small Raa indicate strong interactions with 
the medium and a sizeable energy loss. Experimental results reveal that both quantities are of the 
same order as the respective values for light particles [3-5]. 

The leading order perturbative QCD cross section with a constant coupling as = 0.3 and the 
regular Debye screening mass for the t channel is too small to build up an elhptic flow, which is in 
agreement with the experimental data [14]. However, if we take the running of the coupling into 
account and determine the screening mass from comparison to hard thermal loop calculations, we 
obtain an elliptic flow and Raa, which are much closer to the data. 

The following calculations are done analogously to [15-17]. An effective running coupling is 
obtained from measurements of e^e^ annihilation and non-strange hadronic decays of t leptons 
[15, 18]. Since the t channel of the gQ — )• gQ cross section is divergent, it is screened with a mass 
proportional to the Debye mass mo- 

7 ^ . ^ 2 (4.2) 
t t — Kmf) 

The Debye mass is also calculated with the running coupling for consistency [17]. The prefactor K 



in Eq. (4.2) is mostly set to 1 in the hterature without a sophisticated reason. However, one can fix 
this factor to K" 0.2 by comparing the energy loss per unit length dS/dx of the Born cross section 
with K to the energy loss within the hard thermal loop approach [15-17]. 

This more accurate treatment increases the cross section of elastic gluon heavy quark scattering 
by about a factor of 10. Fig. ^ shows the elliptic flow V2 and nuclear modification factor Raa for 
heavy quarks and for heavy flavor electrons as a function of the transverse momentum pr- To yield 
the same values for these variables as the experimental data, the leading order cross section of 
elastic collisions with the running coupling and improved Debye screening must still be multiplied 
by a K factor of 4. We assume that this artificial K factor stands for the contribution of radiative 
energy loss. However, it must be checked if these corrections have indeed a similar effect as a 
constant K factor of 4. Therefore, the calculation of the next-to-leading order cross section is 
planned for the near future and will complement 2 3 interactions for gluons, which are already 
implemented in BAMPS [10]. 

Especially for high pj, the shape of the V2 curve of heavy quarks is different from the ex- 
perimental data. The reason for this discrepancy is that, experimentally, due to confinement, not 
heavy quarks, but heavy flavor electrons are measured. The latter stem from the decay of D and B 
mesons, which in turn are produced during hadronization of the QGP and which consist of a charm 
or bottom quark and a light quark. However, despite the hadronization and decay processes heavy 
flavor electrons still reveal information about heavy quarks. Essentially, the shape of their spectrum 
is the same as for heavy quarks, but shifted to lower pj due to the decay process. 

For the description of the hadronization process of charm (bottom) quarks to D (B) mesons, we 
use Peterson fragmentation [19]. The decay to heavy flavor electrons is carried out with PYTHIA. 
The theoretical curves for heavy flavor electrons in Fig. ^ are in good agreement with the exper- 
imental data for high pj. For lower pj Peterson fragmentation is not a good description of the 
hadronization and another scheme like coalescence must be employed. In the coalescence picture. 
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Figure 2: Elliptic flow V2 (top) and nuclear modification factor Raa (bottom) of heavy quarks and heavy 
flavor electrons with pseudo-rapidity \ri\ < 0.35 for Au+Au collisions at RHIC with an impact parameter 
of b = 8.2 fm. The cross section of gQ gQ is multiplied with the factor K — 4. For comparison, data of 
heavy flavor electrons [5] is shown. 
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the light quarks of the D/B mesons contribute also to its elliptic flow or nuclear modification factor, 
which increases both. 

As a note, in contrast to previous results [20, 21] we employ here initial heavy quark distribu- 
tions from MC@NLO [22] and not from PYTHIA, since the former employs next-to-leading order 
processes for heavy quark production which describes p-i-p data more accurately [17]. This changes 
the Raa and V2 slightly at intermediate pr. 

Studies on the V2 and Raa of gluons in BAMPS are presented in [23-27]. 



5. // i/A^ suppression 

Two classes of phenomena are important for studying J/y suppression in heavy-ion collisions: 
cold nuclear matter effects and hot nuclear matter effects. The latter are effects that occur due to the 
presence of the quark-gluon plasma which is expected to be produced in ultra-relativistic heavy-ion 
colUsions. Already 30 years ago 7/ i/a suppression due to melting in a medium was predicted to be 
a signature of the QGP [6]. More recently, the regeneration of J/ from a charm and anti-charm 
pair received also much attention and is expected to be important at LHC and maybe also RHIC 
energies [28-34]. 

Cold nuclear matter effects on the other hand are all phenomena ofJ/y suppression that would 
also be present if no QGP was formed. The best systems to study these effects are nucleon-nucleus 
colhsions since one can measure directly the impact of nuclear effects on 7/ i/a production in the 
absence of a QGP. In ultra-relativistic heavy-ion colUsions, however, both effects are present and it 
is a challenge to disentangle them. 

In other words, cold nuclear matter effects influence the initial //i/A production while hot nu- 
clear matter effects describe the modification of the //i/A yield during the evolution of the medium. 

As a note, the presented ingredients of the model and our first results are preliminary. We tried 
to incorporate all important effects in the most accurate way. However, much more work on the 
details has to be done and uncertainties of employed parameters need to be studied. 

The dominant process of initial production is gluon fusion, g + g + Therefore, 

the differential cross section for 7/ i/a production in proton-proton (p-f-p) colhsions is given by 

, J/Xj/ J 

'^^PP - W.(xi,Mf)/,(x2,Mf)^^^^ . (5.1) 



dpTdyj/^dyg dt 

In the present paper we will parametrize the measured J/y production cross section in p-f-p 
colUsions and incorporate cold nuclear matter effects as outUned in the following to get the pro- 
duction cross section for A-i-A colUsions. 

The most important contributions to cold nuclear matter effects are shadowing, nuclear ab- 
sorption and the Cronin effect. Shadowing describes the phenomenon that the parton distribution 
functions of partons in a nucleus f^{x,^F) are modified compared to parton distribution functions 
in a nucleon /™'^^™"(x,jUir). The ratio of both for parton /, 
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can be obtained, for instance, from nucleon-nucleus collisions. In the present study we employ 
the shadowing parametrization EPS08 [35] for Rf{x,flF) and set the factorization scale flp = 
^Pr-^^^lxf, to the transverse mass of the J /\lf. Furthermore, we use a shadowing function that 
also depends on the transverse position in the collision [36] 

^f(xr,x,MF) = l+iVA,p [/?f(x,Mf)-l]^^. (5.3) 

Ta{xt) denotes the nuclear thickness function from the Glauber model andA^Ajp =A Ta{0)/Tab{0) 
is a normalization constant. 

A produced J/y can also be absorbed by the remains of the collided nuclei shortly after its 
production. The nuclear absorption can be effectively described by a survival probability for J/xj/ 
passing through nuclear matter 

= g-<''abs[^i(x7-,ZA,+~)+rA(XT-b,-oo,Zfl)] ^ (5 A) 

Ta {xj- , Zi , Z2) is the nuclear thickness function, but not integrated over the full z range {—°°, +°°) but 
over the interval [zi ,22]- The limits in the formula above are chosen in such a way that it represents 
the path length which the produced J/y travels through the passing remains of the nuclei. For the 
absorption cross section we employ a value of C7abs = 2.8mb [37]. However, the exact value or 
whether nuclear absorption is still present at RHIC and LHC energies is uncertain. 

The J/ Y production cross section in p+p collisions can be parametrized as [38, 39] 

daiC^ _2(„-l)^^^^y"d<^^ 



prdprdy D{y) \ D{y) J dy 

where n = 6 and Z)(j) = {p^)^p{n — 2){\—y^/Y^). 7 = arccosh(^/s^/(2my/i^)) is the the maxi- 
mum rapidity and (p^)pp = 4.14GeV^ the averaged transverse momentum squared. For d<7pp"''/dy 
we employ a double Gaussian distribution [38, 39]. 

To account for the Cronin effect which describes the pr broadening of the fusing gluons in the 
nuclei, we add an additional path length dependence to the mean pj in the parametrization of the 
p+p cross section above: 

{P^) = {P^)pp + agNL (5.6) 

witii agN = 0.1 GeV^/fm [40] and 

L= — [rA(xr, -°°,za) + Ta{xt - h,ZB, +°°)] (5.7) 
no 

which is the path length of the two incoming gluons through the nuclear matter of the other nucleus, 
no denotes the maximum nuclear density from the Woods-Saxon distribution. 
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With all these cold nuclear matter effects the differential 7/ production cross section can be 
written as 



^PT^yj/^dyg 



J dxT J dzA J dzB?i^(xr,ZA)«A(xr-b,ZB)^^(xr,jci,/if)^^(xr-b,jC2,AiF) 



JdXT J dZA J dZBnA{XT,ZA)nA{xT-h,ZB)^fi^T,Xl,HF)^gi^T-h,X2,IJ-F) 

(5.8) 



,-o-abs[rA(xr,^A,+")+rA(xr-b,--zg)] ""PP 



<iPTdyj/^dyg 



where Eq. ^5.1^ was used. In the following we assume that the emitted gluon in the J/y production 
process is soft. As a consequence, xi and X2 are independent of its rapidity yg and we can integrate 
it out: 



^Prdyj/^ 



J dXT J dZA J dZBnA{xT,ZA)nA(xT-h,ZB)^f{liT,Xl,IJ-F)^g{'!iT-h,X2,IJ^F) 

g-CTate[7:i(xr,ZA,+~)+rA(xr-b-=>o,zg)] ""PP ^5 

<ipTdyj/^ 

If there are no cold nuclear matter effects (7?^ = 1, agN = and aabs = Omb) this formula 
simplifies to the well known binary scaling from p+p to A+A collisions: 

^ = TAA{h)p^ (5.10) 
dprdy dprdy 



The J/y distribution obtained with cold nuclear matter effects according to Eq. '5.9' is used 
as an input for our partonic transport model BAMPS, which we employ to study hot nuclear matter 
effects. In the medium, J/y can dissociate via the process //i/A + g— s-c + c whose cross section 
has been calculated to [41, 42] 

/ \3/2 

where w = Pj/^Pgfi/Mj/^, = {s — Mj^^)/2Mj/^ is the gluon energy in the rest frame of the J/y 
and Sj/xf/ = 2Mo —Mjj^ is the binding energy of the J /Y- The back reaction of this process, 
c + c ^ J /Y + g, that is, the regeneration of J/y via charm anti-charm annihilation, is also taken 
into account in BAMPS. The cross section can be obtained from Eq. (5.11^ via detailed balance 

Lattice results [8, 9] indicate that a. J/y can survive in the QGP up to the dissociation temper- 
ature T^i and melts at higher temperatures. In this study we use T^i = 2Tc with = 165 MeV being 
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Figure 3: Raa of J/ y/ at mid-rapidity \y\ < 0.35 in central Au+Au collisions at RHIC as a function of time. 
For comparison experimental data [43] is also shown. 



the phase transition temperature. In BAMPS we implement the effect in the following way. If the 
temperature in a cell is larger than the dissociation temperature T^, all the in this cell decay 
to charm and anti-charm quarks. This is a rather crude and somehow artificial treatment of this 
phenomenon. However, it could be improved by considering, instead of Eq. (5.11)| , a more sophis- 



ticated cross section for J/y dissociation, which leads to such a melting above T^i by itself without 
the need of an additional cutoff. A first attempt in this direction has been done, for instance, in 
Ref. [40] with quasi-free scattering. 

Figure ^ shows the nuclear modification factor Raa of //i// for central Au-i-Au collisions at 
RHIC as a function of time. The initial value is already smaller than 1 which is a consequence of 
cold nuclear matter effects. Then right from the beginning, a lot of initially produced charm quarks 
melt due to the large temperature of the QGP at RHIC. A counter-effect is the J/xj/ regeneration 
which enhances the total J/xj/ number. Our final value is a bit smaller than the experimental data 
point. However, there are significant uncertainties in the initial suppression. In addition, a different, 
more sophisticated regeneration cross section could increase the regeneration and lead to a better 
agreement with the data point. 

In Fig. ^ the Raa of 7/ i/a at RHIC is depicted as a function of the number of participants. As we 
saw in Fig. ^ for central collisions, our final Raa values lie also for non-central collisions slightly 
below the experimental data points, but reproduce the overall shape of the data. This systematic 
underestimation indicates that either our suppression due to cold nuclear matter effects is too strong 
or the regeneration cross section is too small. 
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Figure 4: Raa of y/i// at mid-rapidity \y\ < 0.35 for AuH-Au collisions at RHIC as a function of the number 
of participants, together with experimental data [43]. 



6. Conclusions 

The production and evolution of heavy flavor particles have been studied within the partonic 
transport model BAMPS. Charm and bottom production in the medium produced at RHIC are 
negligible compared to the initial heavy quark yield. At the LHC, however, secondary charm 
production can reach values comparable to the initial yield while bottom production in the QGP 
can also be neglected here. 

The leading order cross section for heavy quark scatterings with particles from the medium 
is too small to explain the experimentally measured elliptic flow and nuclear modification factor. 
However, a more precise Debye screening and the explicit running of the coupling enhances the 
cross section and yields results for heavy flavor electrons, which are much closer to the data, al- 
though a K factor of 4 must be employed for a good agreement with the data. In the future we 
will study if this simple multiplication of the cross section can indeed account for higher order 
contributions. 

Furthermore, we investigated J/y suppression at RHIC. To estimate the initial J/y distribu- 
tion we parametrized the p-i-p production cross section and took for the scaling to A-i-A collisions 
cold nuclear matter effects such as shadowing, nuclear absorption and the Cronin effect into ac- 
count. The space-time evolution of J/y was carried out with BAMPS, which also allows dissoci- 
ation and regeneration of J/y- Preliminary results of the nuclear modification factor Raa of J/y 
obtained with BAMPS for forward and mid-rapidity are slightly smaller than the measured data, 
but resemble the overall shape. Reasons for the smaller yield could be an overestimation of sup- 
pression due to cold nuclear matter effects or the small regeneration cross section. In a future study 
we will investigate this further and also perform calculations for the LHC. 
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